A0'-A068  878  COMPUTER  APPLICATIONS  INC  NEW  YORK  SYSTEMS  ENSINEERI— ETC  F/6  20/1 
INFINITE  ARRAYS  WITH  N0NRI6I0  INTERSTICES. (U) 

SEP  69  V MAN6ULIS  N00024-69>C»1342 


UNCLASSIFIED  CAI-SED-7105-I-Tn-I  NL 


|o^/ 

fiS 68878 

m 

8 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

END 

DATE 

FILMED 

6-79 

ooc 

Incorporated 


miumm  d) 


mart  Program,  Phase.!,  Radiation  Theory 
Technical  Note 


INFINITE  ARRAYS  WITH  NONRIGID  INTERSTICES 


V.  MANGULIS 


CAl/SED-7105  Report  No.  I-TN-1 
Contract  No.  N00024-69-C-1342 


DISnUBUnON  statement  a 


'Nppioved  fo(  pubUo  nleoMi 
Distiibutioii  Unlimited 


E^.'MOST  Project 


MART  Program,  Phase  I,  Radiation  Theory 


Technical jNotc 


JP'INITE  ^RAYS  WITH  JpNRIGID  JNTERSTICES  . 


CAl/SED-7105  Report 
Contract^-MoJ 


Submitted  to 


CODE  NSHP-OOVID 
NAVAL  SHIP  SYSTEMS  COMMAND 
Washington,  D.  C.  20360 


COMPUTER  APPLICATIONS  INCORPORATE: 

V 

Systems  Engineering  Division 


PlSTRaunON  STATEMENT  A 

Approved  foe  public  lelease; 
IHitiibutiait  Unlimited 


ABSTRACT 


Th£  radiation  from  an  infinite  phased  planar  array  of  rectangular  pistons 
(all  with  the  same  velocity  magnitude)  in  a nonrigid  baffle  is  evaluated. 

The  nonrigid  interstices  between  the  array  elements  are  approximated  by 
equivalent  undriven  pistons  with  some  internal  Impedance.  When  the  dimen- 
sions of  the  pistons  are  small  compared  to  a wavelength,  simple  approxima- 
tions are  obtained  for  the  radiation  impedance  of  an  array  element  and  the 

loss  in  source  level  due  to  the  effect  of  the  interstices . ^ The  source  level 

■ 2 ^ 

is  multiplied  by  the  factor  1/  | I + (AA/Zj  cos  0)  | , where  AA  is  the  fraction 
of  the  array  area  occupied  by  interstices,  0 is  the  angle  between  the  normal 
to  the  array  and  the  direction  to  which  the  array  is  phased,  and  Zj  is  the 
interstitial  Impedance  per  unit  area  divided  by  pc,  p = density  of  water,  c = 
velocity  of  sound  in  water. 
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I.  INTRODUCTION 


When  radiation  froM  sonar  arrays  is  considered,  usually  one  prescribes  the 
velocity  distribution  of  the  array  elements,  and  one  assumes  that  the  baffle, 
if  there  is  one,  is  rigid.  With  such  assumptions  the  mathematical  problem  is 
frequently  simpler  than  with  some  other  possible  assumptions.  However,  the 
interstitial  areas  between  the  array  elements  often  are  not  rigid;  in  some 
arrays  there  is  no  baffle  at  all,  in  other  arrays  the  baffle  material  nay 
not  be  rigid  enough.  Nonrigid  surfaces  may  have  a profound  effect  upon  the 
propagation  of  sound. ^ ^ Consequently,  a study  of  the  effects  of  nonrigid 
interstices  is  of  considerable  importance  in  the  design  of  sonar  arrays. 

We  would  like  to  obtain  for  the  effects  of  interstices  a mathematical 

solution  which  is  simple  and  yet  accurate.  The  mathei.Tatlcal  modeling  of 

finite  arrays  on  finite  baffles  is  complicated  even  if  the  baffles  are  rigid. 

However,  the  radiation  of  ^ound  from  an  infinite  phased  array  on  a rigid  plane 

baffle  can  be  described  by  relatively  simple  approximate  expressions  under 

certain  conditions,  and  the  nearfield  sound  pressure  has  been  shown  to  be 

roughly  the  same  for  infinite  and  large  finite  arrays  for  angles  not  too 
7 8 

near  to  endfire.  * . Consequently,  we  have  evaluated  the  effects  of  non- 

rigid  Interstices  in  infinite  arrays;  such  effects  are  expected  to  approxi- 
mate the  effects  in  large  finite  arrays. 

The  nonrigidity  of  the  interstitial  baffle  areas  can  be  modeled  in 
various  ways;  the  simplest  mathe:natlcal  condition  is  to  set  the  pressure 
on  the  Interstice  equal  to  the  velocity  of  the  interstitial  baffle  surface 
times  an  interstitial  impedance  per  unit  area.  Such  a condition  implies 
that  the  nonrigid  interstitial  surface  reacts  only  to  the  local  pressure, 
and  the  movement  of  the  surface  is  not  affected  by  the  acoustic  field  else- 
where. For  example,  if  the  Interstitial  volume  between  the  sonar  array 
elements  is  filled  with  water,  such  cavities  may  have  a wide  range  of 


impedances,  depending  on  the  geometry  of  the  cavities;  in  general  the  cavi- 
ties may  be  coupled  to  each  other,  and  therefore  the  movement  of  the  inter- 
stitial* surface  may  depend  not  only  on  the  local  pressure,  but  also  on  the 
pressure  at  the  other  cavities.  We  will  assume  that  such  coupling  of  cavi- 
ties is  weak  enough  so  that  at  the  interstices  we  can  use  the  local  boundary 
condition:  pressure  equals  velocity  times  an  Impedance  per  unit  area. 
Furthermore,  if  we  break  up  the  interstitial  area  into  smaller  areas  with 
dimensions  much  less  than  a wavelength,  then  each  soialler  area  will  move 
with  roughly  the  same  velocity,  and  therefore  can  be  regarded  as  a piston; 
i.e.,  the  interstitial  area  between  the  active  radiating  pistons  can  be 
broken  up  into  small  equivalent  passive  undriven  pistons  with  some  Internal 
Impedance.  The  radiating  pistons  have  a prescribed  velocity  distribution 
which  is  obtained  by  applying  the  appropriate  voltages  to  the  transducer 
elements.  The  equivalent  interstitial  pistons  are  vibrated  by  the  acoustic 
field  only;  their  velocity  distribution  depends  on  the  mutual  coupling 
between  the  interstitial  and  radiating  pistons,  and  on  the  interstitial 
impedance.  A similar  approach  has  been  used  to  model  a nonrigid  baffle 

9 

surrounding  a small  array  (not  the  interstitial  baffle  between  the  elements) . ’ 

Thus  we  will  consider  an  infinite  array  of  rectangular  pistons  with 
the  interstices  approximated  by  another  set  of  rectangular  pistons  as  shown 
in  Figure  1.  The  radiating  pistons  are  unshaded  in  Figure  1,  and  the  inter- 
stitial pistons  are  shown  shaded. 

All  interstitial  pistons  have  the  same  Internal  impedance  per  unit  area. 

To  be  specific,  we  are  showing  two  interstitial  pistons  for  each  radiating 
piston  in  Figure  1;  however,  if  desired,  the  Interstitial  area  could  be 
divided  further  into  still  simaller  equivalent  pistons. 

In  the  next  sections  we  will  solve  for  the  velocities  of  the  interstitial 
pistons  when  the  velocities  of  the  radiating  pistons  are  specified,  and  we 


will  obtain  the  nearfleld  and  farfleld  pressure  and  the  radiation  Impedance 
of  an  element.  Simple  approximate  expressions  for  the  radiation  impedance 
and  the  loss  in  source  level  will  be  obtained  for  the  case  when  the  dimen- 
sions of  the  radiating  and  Interstitial  pistons  are  much  less  than  a wave- 
length. Numerical  results  will  be  presented  for  the  approximate  loss  in 


source  level*. 


II.  EXACT  EQUATIONS 


Let  the  infinite  array  and  the  nonrigid  baffle  occupy  the  plane  z=0 
(see  Fig.  1).  The  centers  of  the  radiating  pistons  are  located  at  x=ndjj, 
y*indy,  n=o.  ^ I,  ^2,  ...,  m=0,  i I,  ^ 2,  ...  . The  velocity  (in  the  z 
direction)  of  the  radiating  piston  in  the  row  and  the  nf^  column  is 
forced  to  be 

where  vjj  is  a complex  constant,  u)  the  radian  frequency,  t the  time,  and 
=-mkdy  sln9  sin^  - nkdj^  sin0  cos^,  (2) 

where  k=u)/c,  c the  velocity  of  sound  in  water,  and  6,  are  angles  in 
spherical  coordinates  which  determine  the  direction  to  which  the  farfield 
pressure  maximum  is  steered. 

Let  the  velocity  of  the  horizontally  oriented  interstitial  piston 
above  the  ntt  radiating  piston  in  Fig.  I (with  center  at  x=ndx, 

y*(m  + ^)dy)  be  denoted  by 

vJJn  = vjje  itt)t+iYmn,  (3) 


and  let  the  velocity  of  the  vertically  oriented  interstitial  piston  to  the 
right  of  the  n^^  radiating  piston  in  Fig.  1 (with  center  at  x*(n+%)djj, 
y*rody)  be  denoted  by 

For  simplicity  we  will  refer  to  pistons  with  velocities  v&n»  vHn»  or  vXn  as 


1 


R-type,  H-type,  or  V-type  pistons. 


. While  Is  known,  and  are  not,  and  they  have  to  be  obtained 
from  two  equations  for  a force  balance  as  shown  in  Appendix  A.  The 
solution  of  the  two  simultaneous  equations  gives 
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where pcZj  is  the  internal  impedance  per  unit  area  of  both  types  (hori- 
zontal or  vertical)  of  the  equivalent  interstitial  pistons  and  Z^  is 
dimensionless;  p is  the  density  of  water;  and  for  J=R,  H,  or  V,  and  K=R, 
H,  or  V,  the  dimensionless  quantities  S are  given  by 


iAj  D(Aj)  D(Aj^)e^Z(J>K) 

A m = -<»  n = -oo  T(n,  m) 


(7) 


where  A is  the  area  of  the  piston  of  type  J(=R,  H,  or  V),  A = d d = A 

^ c X y ] 

+ A^»  D(Aj)  and  D(Aj^)  are  related  to  the  directivity  of  a rectangular 
piston  of  type  J and  K, 


■><V 


sin  |jg(m)kay]|  sin  [['%h(n)ka^^ 
%g(m)kay  %h(n)ka;J 


(8) 


L 


are  the  dimensions  of  the  J — piston  in  the  x and  y 


direction 


+ 2nn/kd 


g(m)  = sin6  sin^  + 2nm/kd 


E(H,V)  = E(V,H)  = %h(n)kd  - %g(m)kd, 


For  example,  if  we  write  out  explicitly  the  equations  for  J = R, 


K = V,  we  obtain  from  Eq.(7) 


i%h  (n)  kd. 


[%g(m)kSy] 


% g(m)ks. 


[ %h(i:^s  1 sin  f %g(m)kd  1 sin  £ %h(n)kt  1 
X • y * ^ 


%h(n)kt 


%g(m)kd 


%h  (n)  ks 


As  shown  in  Appendix  A,  the  total  pressure  pcv  e P can  be 

R T 


evaluated  as 


P =■  P -.+  (v  /v  )P  + (v  /v  )P 
T R ' H H ' V R-'  V 


where  P^,  Pj^,  P^,  and  P^  are  dimensionless  quantities,  and  for  J = R, 


H,  or  V 


oO  ^ ^ -ih(n)kx  - ig(m)ky  - T(n,m)kz  + iE(J,R) 


where  D(A  ) is  given  by  Eq.(8),  and  E(J,R)  is  given  by  Eq.(13). 

U 

The  radiation  impedance  of  a radiating  piston  is  pcs  s Z , where 

X y R 

the  dimensionless  quantity  Z is  derived  in  Appendix  A, 

R 


Z = S + (v  /v  )S  + (v  /v  )S  , 
R RR  H R HR  V R'^  VR’ 


and  and  are  given  by  Eq.(7). 


III.  APPROXIMATIONS 

While  the  exact  equation  for  the  pressure  field  contains  an 
infinite  number  of  terms  as  in  Eq.(16),  it  has  been  shown  that  for 
small  pistons  close  together  a single  term  (the  m=0,  n=0  term)  is 
sufficient  to  approximate  the  exact  equation.  Thus,  if  the  array 
spacings  d and  d are  much  smaller  than  a wavelength,  then,  as  shown 

X y 

in  Appendix  B,  we  obtain  the  following  approximate  expression  for 
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f 

ii: 


y 

[i 


the  total  radiated  pressure 


P « G.  R jexp  (-Ikxslndcos^  - IkyslnOsln^ 
T A^cosQ 


• -ikzcosQ)  = G.P^, 

and  the  radiation  impedance  of  an  element  becomes 


(18) 


Z « , 

R A cos6 


G ~ • G » 


(19) 


where  P^  and  Z^  are  the  approximate  pressure  and  impedance  which 
one  would  have  obtained  for  small  pistons  close  together  if  the 
interstices  had  been  rigid. ^ Consequently,  G is  the  factor  which 
gives  the  effect  of  the  nonrigid  interstices. 


Z^cosQ 


Z^cosS /AA 


AA  + Z^cosQ  1 + (AA/Zj.cos0)  1 + (Zj.cos6/AA) 

(20) 

where  AA  is  the  fraction  of  the  array  area  occupied  by  the  inter- 
stices, 


AA  - (A  + A ) / A 
H V C 


(21) 


Note  that  G approaches  1 as  it  should  when  either  AA  -•  0 or  |Zj.  | -♦  *. 

For  AA  “•  1 G is  finite,  but  P and  Z vanish  because  A -•  0.  G “•  0 

r r R 

if  Zj  “•  0,  i.e.,  the  array  does  not  radiate  if  the  interstices  are  per- 

fect pressure  release  surfaces  (at  least  in  this  approximation  the  pressure 
vanishes;  there  may  be  some  of  the  smaller  terms  in  the  exact  equation 
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[; 


t 

I 

I 


I -- 


i I 


0 


which  do  not  vanish) . Observe  that  G depends  only  on  the  combination 
ZjCogO/AA,  not  on  Z^.,  6,  or  AA  separately.  Z^  and  AA  are  dimension- 
less, therefore  Z^cos0/AA  is  also  dimensionless. 

Vfhen  8 -*  90®,  then  G ”•  0,  but  -*  » and  Z^  -♦  » and  the 

products  GP  and  GZ  remain  finite.  However,  the  behavior  of  the 
r r 

infinite  array  does,  not  approximate  the  behavior  of  the  finite  array 
for  angles  9 too  close  to  90°. 

Note  that 


A„  / A^  cos6 
R C 


(22) 


is  a pure  resistance,  and  under  the  right  conditions  the  presence  of 

the  nonrigid  interstices  will  make  Z complex.  We  are  mainly  inter- 

R 

ested  in  the  change  in  the  power  output,  proportional  to  the  radiation 
resistance,  and  thus  proportional  to  the  real  part  of  G, 


Re  (G) 


where  R^  + iX^ 


2 2 2^ 

(Rj;  + Xj)cos  0 + AARjCosO 

(AA  + Rj.cos6)  + (XjCOsS)" 


(23) 


On  the  other  hand,  the  source  level  in  the  far  field  is  propor- 
tional to  thus  it  is  proportional  to  [cp. 


(R^  + X^)cos^e 


(AA  + RjCosS)^  + (XjCosS) 


(24) 
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The  difference 


Re 


(G)  - I g| 


LA 


RjCOs6 


(A  A + RjCosO)^  + (XjCos®)^ 


(25) 


Is  proportional  to  the  power  dissipated  in  the  interstices;  the  differ- 
ence vanishes  when  = 0 because  when  the  interstice  impedance  is 

purely  reactive  then  the  interstices  cannot  absorb  power. 

1 2 

Figure  2 shows  J Gj  vs.  | Z^cos6/AA  j for  three  different  conditions: 

1)R  =0;2)R  =X;3)X=0.  The  source  level  is  highest  for 

I I I I 

R = 0,  lowest  for  X = 0;  i.e.,  for  the  same  value  of. the  interstice 

I • . . ■ I 

impedance  magnitude,  a pure  resistance  produces  a greater  loss  of  source 
level  than  a pure  reactance  because  the  interstice  then  not  only  modifies 
the  radiated  pressure  but  also  absorbs  power.  Note  also  that  |g|^  is 
independent  of  the  sign  of  X^. 

IV.  CONCLUSIONS 

We  have  derived  the  exact  expressions  for  the  pressure  [Eq.(15)] 
and  radiation  impedance  |j:q.(17)J  for  an  infinite  array  of  rectangular 
pistons  with  passive  interstitial  pistons  interspersed  among  the  radi- 
ating pistons.  For  simplicity  we  assumed  that  the  nonrigld  interstices 
can  be  approximated  by  two  equivalent  pistons  for  each  radiating  piston; 
however,  we  could  have  specified  without  much  difficulty  more  intersti- 
tial pistons  for  each  radiating  piston.  With  two  interstitial  pistons 
we  had  to  solve  two  simultaneous  equations;  with  N interstitial  pistons 
we  would  have  to  solve  N simultaneous  equations. 

Due  to  a lack  of  funds  we  have  not  evaluated  the  exact  equations 
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numerically.  However,  for  small  pistons  close  together  we  have  been 

able  to  obtain  a very  simple  approximate  expression  [^Eq.(20^  for  the 

factor  which  modifies  the  pressure  and  radiation  Impedance  because  of 

the  presence  of  thenonrlgld  interstices.  The  approximate  expression 

shows  tha't  for  angles  near  broadside  the  source  level  will  be  reduced 

by  less  than  4 db  if  the  Interstices  occupy  less  than  50%  of  the  array 
« 

area  and  if  the  magnitude  of  the  normalized  dimensionless  Interstitial 
Impedance  |Zj[  is  greater  than  I (see  Figure  2).  Consequently,  one 
should  be  able  to  design  an  array  with  Interstices  In  such  a way  that 
the  detrimental  effects  of  the  Interstices  are  tolerable  — if  the 
presence  of  pressure  release  materials  and  resonances  In  Interstitial 
cavities  is  avoided,  then  one  should  be  able  to  riiake  jZ^j greater  than  1. 

The  reduction  of  the  area  occupied  by  interstices  would  also  help. 

APPENDIX  A. 

DERIVATION  OF  THE  EXACT  EQUATIONS 

#"h 

since  the  piston  velocities  and  the  pressure  field  at  the  m^  row 

and  the  n—  column  will  be  the  same  as  at  m = 0 and  n = 0,  except  for 

the  phase  Y , we  can  confine  our  attention  to  the  pistons  at  m = 0,  n = 0. 
tnn 

First  consider  the  pressure  field  due  to  the  radiating  pistons  alone, 
l.e.,  assume  the  interstices  are  rigid.  Then  the  pressure  field  can  be 
obtained  from  Eq.(14)  in  Reference  7;  in  Reference  7 the  pistons  are  circu- 
lar, therefore  the  integration  in  Eq.(14)  of  Reference  7 is  over  the  circu- 
lar area  of  a piston;  we  must  replace  that  by  integration  over  the  rectangu- 
lar area  of  our  piston,  which  then  yields  the  dimensionless  quantity  P as 

R 

given  by  Eq.(16)  in  Section  II  here. 


Next  let  us  calculate  the  pressure  pev  e^^P  when  the  H-type 

H H 

pistons  vibrate  with  velocities  v^  and  the  rest  of  the  xy-plane  Is 
• mn 

rigid.  P is  given  by  the  same  equation  as  P , if  we  replace  the 
**  R 

dimensions  of  the  radiating  pistons  by  the  dimensions  of  the  H-type 
pistons,  and  if  we  compensate  for  a shift  of  origin  from  the  R-type 
to  H-type  pistons;  the  latter  accounts  for  the  phase  E(H,R)  in  Eq.(l6). 

P^  is  obtained  the  same  way. 

Thus,  if  we  know  the  velocities  v , v , v , then  the  total  pressure 

R H V 

is  obtained  from  the  superposition  of  the  three  partial  pressures  due  to 
each  type  of  pistons: 

iwt„  iODt  , 

pev^e  Pj.  p«  <Vr  Vh 

(A-l) 

The  superposition  yields  the  correct  expression  because,  for 

example,  P^  gives  the  correct  velocities  for  the  H-type  pistons,  and 

It  adds  nothing  to  the  velocities  of  the  R-type  or  V-type  pistons, 

since  Pjj  is  a solution  of  the  model  where  only  H-type  pistons  vibrate, 

and  R-type  and  V-type  pistons  are  immobile.  Similar  considerations 

apply  to  P and  P . 

R V 

However,  this  far  we  still  have  not  found  the  relationship 

between  v and  v , v . To  obtain  that  relationship,  we  must  consider 
R H V 

the  forces  on  the  H-type  and  V-type  pistons,  i.e.,  the  boundary  condi- 
tion on  the  interstices. 

To  find  the  total  force  on  an  H-type  piston,  we  must  integrate 

the  total  pressure,  given  by  Eq.(A-l),  over  the  area  A . Define 

H 

S . i_  ‘/JdA  P 
JK  Aj^  A^  J, 

<A-2) 
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i.e.,  the  dimensionless  quantity  S is  the  Integral  of  the  partial 

«J1C 

pressure  due  to  all  J-type  pistons  over  the  area  of  one  K-type  piston, 
divided  by  the  area  A ; the  integration  yields  Eq.(7)  In  Section  II. 

Iv 


Sjjj  1®  proportional ■ to  the  force  on  one  K-type  piston  due  to  the  com- 
bined action  of  all  J-type  pistons.  Then  the  total  force  F on  the 

H 


area  A^ 


is 


F„  ■ pcA  ,e^^  (v^S^  + V S + v S ) 

H ■ ^ R RH  H HH  V 

(A-3) 

However,  if  the  internal  impedance  of  the  H-type  piston  is 
PcA^Zj,  then  the  boundary  condition  on  the  Interstice  (force  * velocity 
times  Impedance)  becomes 


H 


- V e . pcA  Z 
H HI 


(A-4) 


where  the  minus  sign  arises  from  the  fact  that  v is  defined  as  the 

H 

velocity  into  the  water,  while  the  boundary  condition  contains  the 
opposite  velocity  (into  the  equivalent  piston).  If  we  equate  Eqs.(A-3) 
and  (A“4),  we  obtain 


+ vS  + vS  “-  vZ. 
R RH  H HR  V VH  H I 


(A- 5) 


Similarly,  from  the  total  force  on  one  V-type  piston  with 

internal  Impedance  pcA  Z we  find 

V I 


'^R®RV  \®HV  ^ 


- Vi 


(A-6) 


If  we  solve  simultaneously  Eqs.(A-5)  and  (A-6),  then  we  obtain 


I 


Eqs.(5)  and  (6)  In  Section  II.  With  and  v specified,  the  total 

pressujre  as  given  by  Eq.(A-l)  Is  now  completely  determined. 

Moreover,  from  the  total  force  on  one  radiating  piston  we  obtain 

the  radiation  Impedance  pcA  Z , because  by  definition 

R R 


F - pcA„Z  .V  e^*^, 
R R R R * 


(A- 7) 


and,  furthermore,  from  the  Integral  of  the  total  pressure  on  the  radi- 
ating piston. 


F - pcA^e^^^^^Cv  S + vS  + vS). 
R ^ R RR  H HR  V VR^’ 


(A- 8) 


consequently.  If  we  equate  Eqs.(A-7)  and  (A-8),  we  obtain  Eq.(17)  for 

the  normalized  dimensionless  radiation  Impedance  Z . 

R 


APPENDIX  B 

DERIVATION  OF  THE  APPROXIMATE  EXPRESSIONS 
If  kd  -*  0 and  kd  -*  0 (and  then  ks  -♦  0,  etc.,  because  s < d , etc.), 
then  from  Eqs.(9)  and  (10)  h(n)  -*  » and  g(m)  -♦  »,  except  for  n = 0 and  m =0. 
Furthermore,  T(n,m)  -*  «,  except  that 


T(0,0)  = lcos8 


(B-1) 


Consequently,  for  small  pistons  close  together  the  dominant  terms  In  Eqs.(7) 
or  (16)  are  the  terms  with  m = 0,  n = 0.  For  m “ 0,  n * 0,  and  small  pistons 
close  together  the  directivity  D(Aj)  w 1 in  Eq.(8),  and  the  E(J,K)  -*  0. 


(B-2) 


(Aj/A^co80)  exp  (-IkxslnQcos^  - iky8ln08inj^  - ik2CO80) 

(B-4) 

If  we  8ub8titute  Eq8.(B-2),  (B-3),  and  (B-4)  into  Eq8.(15)  and  (17), 
then  we  obtain  Eq8.(18)  and  (19). 
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